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Abstract. 

Electron-impact ionization cross sections for the ls2s and ls2s 
metastable states of Li"*" are calculated using both perturbative distorted-wave 
and non-perturbative close-coupling methods. Term-resolved distorted- wave 
calculations are found to be approximately 15% above term-resolved ij-matrix 
with pseudostates calculations. On the other hand, configuration-average time- 
dependent close-coupling calculations are found to be in excellent agreement with 
the configuration-average H-matrix with pseudostates calculations. The non- 
perturbative ij-matrix and close-coupling calculations provide a benchmark for 
experimental studies of electron-impact ionization of metastable states along the 
He isoelectronic sequence. 
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1. Introduction 

Electron-impact ionization of low-charged ions is an important component of many 
processes in astrophysics and controlled fusion plasmas. This topic has therefore 
generated considerable interest, both experimental and theoretical. Previous studies 
of ground-state He (Pindzola & Robicheaux 1998, Pindzola et al. 2004), as well 
as He-like Li (Pindzola et al. 2000) and Be (Colgan et al. 2003), have shown that 
both the i?-matrix with pseudo-states (RMPS) method (Bartschat 1998) and time- 
dependent close-coupling (TDCC) methods (Pindzola & Schultz 1996, Pindzola & 
Robicheaux 1996) give cross-sections that are in very good agreement with experiment. 

Far less work has been done on the electron-impact ionization cross-section of 
metastable ions. One important example is excited He (ls2s ^S)^ that has been 
studied both in experiment and theory. The experiment, performed by Dixon et al. 
(1976), is in reasonable agreement at high energies with Bethe-Born calculations that 
pre-date the experiment (Briggs & Kim 1971), and is in good agreement with first- 
order plane- wave Born results generated shortly afterwards (Ton-That et al. 1977). 
More recently, however, several non-perturbative calculations have shown a large 
discrepancy with experiment (Colgan & Pindzola 2002, Bartschat 2002, Fursa & 
Bray 2003). These papers used TDCC, convergent close-coupling, and i?-matrix with 
pseudostates (RMPS). All of these calculations are in agreement, and all differ from the 
Dixon experiment by a factor of two, even at high energies. The helium discrepancy 
motivates us to see whether the situation improves with metastable He-like Li. 

Further motivation comes from plasma modeling. It has been previously shown 
that ionization from the metastable Li"*" can dominate the effective ionization rate 
coefficient even at relatively low electron densities (10^" cm~'^) (Loch et al. 2004). 
However, all existing data on the excited states of Li"*" is based on distorted-wave 
calculations, which can get progressively worse as the term energy gets higher; this 
was demonstrated in H-like ions, including Li^+, see Griffin et al. (2005). Thus it 
is important to test the accuracy of the ionization data from excited states and to 
determine the effect that such data has on plasma modeling. 

In this paper we build upon a previous study of the Li+ ground state (Pindzola 
et al. 2000) to examine the electron-impact ionization cross-section of metastable 
ls2s and ls2s Li+. We compare the perturbative distorted- wave method with the 
non-perturbative RMPS and TDCC methods. The cross-sections calculated using the 
non-perturbative methods converge to the same value in Li^ , while the distorted- wave 
method gives results that are significantly higher. This study provides a benchmark for 
future experimental studies of electron-impact ionization from metastable states along 
the He isoelectronic sequence. In the following section we discuss the various methods 
we have used to calculate electron- impact ionization cross-sections of metastable Li+ . 
The results of these calculations are presented and compared in Section [3l 

2. Theory 

2.1. Distorted-wave method 

The distorted-wave (DW) cross sections are calculated from a triple partial-wave 
expansion of the first-order perturbation theory scattering amplitude. Both direct 
and exchange components of the scattering amplitude are included. The incoming 
and outgoing electrons are calculated in a potential, while the ejected electron 
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is calculated in a V^~^ potential, where N = 2 \s the number of electrons in the 
target ion (Younger 1985). The DW potential for all continuum electrons is made 
up of a configuration-average Hartree potential for the direct interaction, and a local 
semiclassical approximation for the exchange interaction. 

For the tcrm-rcsolvcd distorted-wave calculations, the bound orbitals are 
calculated using Fischer's multi-configuration Hartree-Fock code (Froese Fischer 
et al. 1997). A double configuration (Is^ and ls2s) calculation for the ls2s term 
yields an energy of -5.036 a.u. and an ionization potential of 14.59 eV, in good 
agreement with the experimental value of 14.86 eV (National Institute of Standards 
& Technology 2006). A single configuration (ls2s) calculation for the ls2s term 
yields an energy of -5.109 a.u. and an ionization potential of 16.58 eV, in very good 
agreement with the experimental value of 16.60 (National Institute of Standards & 
Technology 2006). 

For the configuration- average DW calculations, the bound orbitals are calculated 
using Cowan's (1981) Hartree-Fock code. Calculation for the ls2s configuration yields 
an energy of 16.07 eV, in good agreement with the experimental value of 16.16 eV. 



2.2. R-matrix with pseudostates method 

The codes used in this work are based on the serial codes published by Bcrrington 
et al. (1995), with modifications for the addition of pseudostates (see, e.g. Gorczyca 
& Badnell 1997), and parallelization described in Mitnik et al. (1999), Mitnik 
et al. (2003), and Ballance & Griffin (2004). In this method, the high Rydberg 
states and the target continuum are represented by a set of orthogonalized Laguerre 
radial wavefunctions: the so-called pseudostates. The ionization cross-sections are 
determined by summing over excitations above the ionization threshold, including all 
single-electron excitations to the pseudostates as well as doubly excited states. 

Our RMPS basis used spectroscopic orbitals up to n = 3 and pseudostates from 
n = 4 to n = 14, with a maximum angular momentum of Z = 4. The pseudostates with 
I = 4 and 11 < n < 14 were omitted to keep the size of the calculation manageable. 
For the incoming electron, partial waves from Z = to 13 were calculated using the 
i?-matrix method with exchange included. This was topped up from I = 14 with 
methods described by Burgess ct al. (1970) (see also Burgess 1974). 



2.3. Time- dependent close-coupling method 

We use the "frozen core" approximation that was previously used in electron-ionization 
from metastable helium, where we freeze the Is electron (see Colgan & Pindzola (2002) 
and references therein). The frozen Is orbital is the ground state of the hydrogenic 
Li^+ ion. A complete set of orbitals are then obtained by diagonalization of the single- 
particle Hamiltonian, including the direct term of the Hartree potential, and a local 
approximation to the exchange interaction. The local exchange potential was scaled 
to give single-particle energies close to experiment. 

The initial two-electron wavefunction, P^^^'^ (ri , r2 , i = 0), is an antisymmetrized 
product of a radial wavepacket and the target 2s radial orbital, with a particular LS 
symmetry. The propagation in time is governed by the Schrodinger equation, which 
we write as a set of time-dependent close-coupled partial differential equations 

^^^%^=Tw.(n,r.)P^£(n,r„t) 
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Figure 1. Term-resolved ls2s and ls2s electron-impact ionization cross- 
sections, ""ionization' S'S^'iiist incident electron energy, E. We show the raw RMPS 
output (light solid lines) along with fits from Equation ((2j (heavy lines). The raw 
distorted-wave data is also shown (circles). (1Mb = lO""*^* cm^) 

where Ti-^i^{ri,r2) is the single particle Hamiltonian, which includes kinetic energy, 
nuclear, direct and local exchange operators, while J/;^;^ ;/ ;/ (ri, couples the 
(hh) scattering channels. Some time after the collision, the two-electron radial 
wavefunctions are projected onto products of the single-particle orbitals to calculate 
the probability of excitation. The ionization probabilities, and hence the ionization 
cross-sections, are found by subtracting the probability of any electrons being bound 
from unity. 

3. Results 

The term resolved DW and RMPS results are shown in Fig. [1] We see that for both 
methods the cross-sections are larger for the ls2s than the ls2s across the entire 
energy range. Furthermore, the DW results are approximately 15% higher than the 
fitted RMPS results. The data fits were made using the formula of Younger (1981): 

f^ionization = ^^(^(l " V«) + ^(1 - 

+ Cln(u) +Dln(u)/u) (2) 

where / is the ionization energy, E is the incident electron energy, and u = E/I. The 
coefficients A, B, C, and D are determined from a least-squares fit to the calculated 
cross-section (fitting parameters are available from the authors upon request). Note 
that the coefficient C can be independently determined from the photoionization cross- 
section, but we have left it as a free parameter. 

While the DW and RMPS calculations are term resolved, and thus return both 
ls2s and ls2s cross sections, the TDCC calculation is configuration averaged 
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Figure 2. Configuration-averaged electron-impact ionization cross-sections, 
"^ionization' against incident electron energy, E. The distorted- wave calculation 
(dashed line) gives a cross-section around 15% larger than the non-perturbative 
RMPS calculation (solid line). Three TDCC points are also shown (squares), and 
these are in good agreement with the RMPS calculation. (1Mb = 10^^® cm^) 

and only gives a single cross-section for the ls2s configuration. In order to compare 
the various theoretical calculations we show all of the configuration-averaged results in 
Fig. [21 The RMPS and results have been converted to a configuration-average 
cross section by averaging the fits of ([2]). The DW cross-section in this graph was 
calculated using the configuration-average approximation (see Section 12.11) . We also 
took an average of the term resolved DW calculations, and it is in good agreement 
with the configuration-average result. The TDCC and RMPS cross sections are in 
excellent agreement with each other, while the DW result is around 15% higher at the 
peak than the non-perturbative results. 

It is interesting that the perturbative DW method produces cross-sections that are 
in good agreement with TDCC and RMPS for the ground state (Is^ ^S") ionization 
(Pindzola et al. 2000) and overestimates the metastable cross section. This is not 
unexpected for such a near neutral system, since DW doesn't include higher-order 
coupling between the outgoing electrons. This "three body" effect is more important 
when the ionization limit is lower, as it is in the case of metastable Li"*". This trend 
is similar to that found by Griffin et al. (2005) for hydrogenic systems, including 
Li^"*". The discrepancy between the perturbative distorted- wave method and the non- 
perturbative RMPS and TDCC methods is therefore larger for ionization from excited 
states than for ground states. 

One would also expect the effect to be more important when the effective ion 
charge (the charge that the escaping electrons see) is smaller. In fact, this trend can 
be seen when one compares the metastable He ionization cross-sections of Colgan 
& Pindzola (2002) with the He-like Li results of this paper. As mentioned in the 
introduction, the DW and perturbative methods differ by a factor of two for metastable 
He. Furthermore, this trend can be observed by comparing ionization from H to 
ionization from Li^+ (Griffin et al. 2005). 
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4. Conclusion 

We present results of three different calculations of the electron-impact cross section of 
metastable ls2s and ls2s Li+: distorted- wave, ii-matrix with pseudostates, and 
time-dependent close-coupling. We find that the non-perturbative methods, i?-matrix 
and time-dependent close-coupling, are in excellent agreement. The perturbative 
distorted- wave method, however, gives cross-sections that are significantly larger than 
those given by the non-pcrturbativc methods. The non-perturbative i?-matrix and 
close-coupling calculations provide a benchmark for future experimental determination 
of absolute ionization cross-sections from metastable states along the He isoelectronic 
sequence. 

The fact that the DW and TDCC methods are actually in agreement for ionization 
cross-section from the Li+ ground state (Pindzola et al. 2000) has previously been 
taken as a sign that DW is fairly dependable for this species. However higher- 
order coupling between the outgoing electrons is not taken into account in the DW 
calculations. Although these correlations become smaller as the nuclear charge Z 
increases, they become larger for excited-state ionization. Thus the non-perturbative 
methods are more appropriate for excited-state ionization. All existing electron-impact 
ionization data for Li+ excited states comes from distorted-wave calculations (Loch 
et al. 2006). In light of the results presented in this paper, however, this data will 
need significant revision. 
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